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The magnetic flux threading a conventional superconducting ring is typically quantized in units 
of $o = hc/2e. The factor 2 in the denominator of <3>o originates from the existence of two different 
types of pairing states with minima of the free energy at even and odd multiples of <3>o- Here we 
show that spatially modulated pairing states exist with energy minima at fractional flux values, in 
particular at multiples of $o/2. In such states condensates with different center-of-mass momenta of 
the Cooper pairs coexist. The proposed mechanism for fractional flux quantization is discussed in the 
context of cuprate superconductors, where /ic/4e flux periodicities as well as uniaxially modulated 
superconducting states were observed. 
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In superconducting (SC) rings flux quantization in 
units of $o = hc/2e is usually attributed to the charge 2e 
of Cooper pairs carrying the supercurrent. This connec- 
tion has been anticipated by F. London [1] long before 
the experimental confirmation of the SC flux quantum $o 
by Doll and Nabauer [2] and by Deaver and Fairbank [3] 
in 1961. In the same year, $ was derived from the BCS 
pairing theory as a consequence of the existence of two 
classes of SC wave functions with energy minima at even 
or odd multiples of i>o [4-6]. 

The formation of a pair condensate alone does not in- 
evitably imply that $o is the unit of flux quantization. 
In fact, this conclusion is only valid for a uniform con- 
densate of non-interacting Cooper pairs, as assumed in 
BCS theory. Normal persistent currents can as well sus- 
tain a hc/2e flux periodicity in special geometries [7] 
or due to electron-electron interactions [8, 9]. Correla- 
tions between Cooper pairs may also lead to fractional 
flux quanta. Such effects were first discussed by Little 
in 1964 in the context of fractional flux periodicities in 
the critical temperature of conventional SC cylinders [10]. 
More recently, unusual flux periodicities were recently re- 
ported in SQUID experiments with high-T c cuprates [11- 
13]. An example for &o/2 oscillations of the critical 
current is shown Fig. 1. The low magnetic field data 
indicate a sin2cp Josephson relation. As one likely ori- 
gin of the fractional periodicity, multiple Andreev scat- 
tering in grain-boundary Josephson junctions was pro- 
posed [11, 12]; however, the abrupt disappearance of the 
$o/2 periodicity beyond a threshold field has remained 
unexplained. In the search for an alternative origin, in- 
teractions between Cooper pairs and quartet formation 
have been investigated theoretically [14-17]. Yet, their 
influence on flux quantization is unresolved. 

Here we formulate a conceptually distinct mechanism 
for fractional flux quanta in superconductors with a spa- 
tially modulated SC order parameter (OP), in particular 



for $0/2 flux periodicity and a sm2ip Josephson relation 
in SQUIDs. The proposed concept rests on the coexis- 
tence of pair condensates with different center-of-mass 
momenta (COMM) of Cooper pairs and is unique for su- 
perconductors with unconventional pairing symmetries. 

Consider the uniform OP A(0) = Ae^-*/* ) for 
electron pairs with quantized COMM hq and angular co- 
ordinate on a hollow cylinder threaded by the magnetic 
flux <I>. In this geometry, q refers to the phase winding 
number of the OP upon circulating the cylinder once. 
The free energy F of this system has a minimum for a 
vanishing phase gradient, i.e., at the flux value $ = <7$o- 
Since q is an integer, the fluxoid threading the cylinder 
is quantized in units of $0 [18, 19]. 

If instead the magnitude A of the OP oscillates around 
the cylinder, Cooper pairs with different COMMs exist. 
An example for a sign changing OP is 

A(0) = (A ie 4910 + A 2 e iq2 °) e"^*/* 

2Ai cos (%! - q 2 }/2) e «([«i+»]/2-»/*o). (1) 



60 




-20 -10 10 20 

poH ([iT) 

FIG. 1: Critical current J C (H) through a 24° YBCO grain- 
boundary SQUID at T = 4.2 K from Ref. [12]. Clearly visible 
is the abrupt change of periodicity at fioH « ±5 [iT. For the 
SQUID size in this experiment, a flux of $0 is achieved for a 
magnetic field of 2.7 uT. (Courtesy of C. Schneider) 
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In such a state, pairs with COMM qi coexist with pairs 
with q 2 ; its phase-winding number is [qi + 92] /2 which 
implies that F has a minimum at <&/$o = [<7i + 92] /2. 
Since q\ + q 2 can be both even or odd, the flux is hence 
quantized in units of $o/2. 

The above state generalizes the time-reversal sym- 
metric "pair-density wave" (PDW) concept with qi = 
—qi (see e.g. Refs. [20, 21] and references therein). 
This state bares similarities to Fulde-Ferrell-Larkin- 
Ovchinnikov states in an external magnetic field [22, 23]. 
Due to its concomitant charge-density modulations, the 
PDW was proposed as a candidate groundstate of striped 
cuprate superconductors [21, 24]. For unconventional 
pairing symmetries with gap nodes [25] it was indeed 
verified that the PDW state has a finite range of sta- 
bility. A Ginzburg-Landau analysis furthermore demon- 
strated that the PDW state can sustain a 3>o/2 vortex 
phase [20, 26]. The melting of the PDW state may 
even give rise to a state with charge-4e superconductiv- 
ity, which also leads to a $o/2 periodicity in a SQUID 
geometry [27]. 

Another system where the generalized PDW concept 
may lead to $o/2 flux quanta is the spin-triplet super- 
conductor Sr 2 Ru04. Here it is the two-component order 
parameter {A|J,A^} for the s z = 1 and the s z = —1 
condensates which allows for the coexistence of differ- 
ent COMMs [28]. Indeed, experimental evidence was re- 
ported that the flux through a Sr2Ru04 ring is quantized 
in units of $0/2 [29]. The rare observations of fractional 
flux quanta have been encountered only for unconven- 
tional superconductors. This is not an accident, because 
it is the existence of gap nodes, in particular for the d- 
wave pairing symmetry of the cuprates, which allows for 
the emergence of a PDW groundstate [25] . 

In the following we apply the concept of the generalized 
PDW with arbitrary combinations of COMMs q\ and q 2 
to flux-threaded hollow cylinders. We show that for suffi- 
ciently large cylinders, the groundstate energy E(n^o/2) 
(with integer n) is minimized by a combination of two 
COMMs qx and q 2 fulfilling (gi + q 2 )/2 = n/2, i.e., the 
flux is quantized in units of $0/2. The energy minima 
for different n become degenerate in the thermodynamic 
limit, leading to a $o/2-periodicity of the groundstate 
energy and the supercurrent. 

We start from a pairing Hamiltonian on a square lattice 
(the lattice constant is set to 1) with periodic boundary 
conditions, and Ng sites in ^-direction, N z sites in z- 
direction [41]: 



ft = X/ k( ^) C ks C k S + 
k,s 

EEE y ( k ' k '' q) C L C -k+q 5 ' C -k'+qs'Ck' S (2) 
q k,k' s.s' 

with an attractive interaction V(k, k',q). Here, k = 
(kg,k z ), where kg — 1,...,N$ enumerates the angu- 



lar momenta Hkg for motion around the cylinder, and 
k z = 1,...,N Z the momenta hk z in z-direction. The 
single-electron dispersion takes the form 



e k (y>) = -2t 
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with nearest-neighbor hopping amplitude t. The mag- 
netic flux enters in e^ip) through ip = &e/hc = $/2$ - 
For a superconducting state with singlet pair- 
ing, we use the BCS mean-field decoupling 
scheme and approximate (c k1 .ct k+q |C_k'+q4.Ck'-|-) — > 

( c kt c -k+qi)( c -k'+qi c k't)- Tne Heisenberg equations of 
motion for the spin independent imaginary time Green's 
function £(k,k',r) = -(T T c ks (T)cj c , s (0)) and the 
anomalous propagators F(k, k',r) = (T T Ck s (r)c_k's' (0)) 
and J"*(k,k',r) = (T t c^ s {t)cI, s , (0)) for s ^ s' then 
lead to the Gor'kov equations [30]: 



£/(k,k',u; n ) = £/o(k,w n ) 



<W-^A(k,q).F*(k-q,k>„) 
q 



(4) 



and 



F(k, k',w n ) = G (k, uj n )J2 A(k, q}g{-k', k + q, -u n ), 

q 

(5) 

where Qo(k, oj n ) = [iu> n — £k + n] 1 1S the Green's func- 
tion in the normal state and co n = (2n — X^-nk^T is the 
fermionic Matsubara frequency. The average charge den- 
sity p is controlled by the chemical potential /1. The 
order parameter A(k, q) represents electron pairs with 
center-of-mass momentum fiq and is determined self con- 
sistently as shown below. 

At this point the form of the interaction is crucial. 
As the simplest ansatz that allows for unconventional 
pairing, we choose an attraction between electrons on 
neighboring sites. With the restriction to singlet pairing 
only the extended s- and the d-wave channels remain, 
which leads to the interaction V(k, k', q) = V+(k, k', q) + 
y_(k,k',q) with 

V ± (k,k',d) = Vg ± (k - q/2) 5± (k' - q/2), (6) 

where g±(k) = cos(2nkg/Ng) ± cos(27rfc z /A z ). The OP 
thus becomes 

A(k, q) = A s (q). 9+ (k - q/2) + A d (q) ff _(k - q/2) (7) 

with a non-vanishing extended s-wave contribution 
A s (q) if q 7^ [25]. Cooper pairs moving in the flux 
threaded cylinder acquire angular momentum and we 
therefore choose q = (q, 0) along the ^-direction. In- 
serting V(k, k',q) from (6) into (4) generates a set of 
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FIG. 2: Energies E q - q+1 (l/4)/(N e N z ) (blue), E q - q (0)/(N g N z ) (pink), and E q - q (l/4)/(N e N z ) (grey) per lattice site as a 
function of the circumference No of the cylinder, for N z — 80, V\ = 2.2 t, p = 0.8 and t' = — 0.3t. Although E qitq2 (tp) oscillates 
with the cylinder's radius, for large enough cylinders E q ,- q (0) = £ 9 ,_ g +i(l/4) is always smaller than E qt - q (l/4). 



coupled self-consistency equations for A s (q) and Ad(q): 

A S)rf (q) = - ^ X! SM(k - q/2) J] ^ - <L k, Wn ). 

(8) 

In an ansatz for a self-consistent solution, we choose Q 
trial vectors qi, . . . , q<g and set A(k, q) = for all q ^= q^. 
Thereby we test selected combinations of q- vectors. This 
procedure and its solutions without magnetic flux are dis- 
cussed in Ref. [25]. In the context of flux quantization, 
two notions are important: (i) If instead of Eq. (6) con- 
ventional s-wave pairing originating from an on-site in- 
teraction is considered, the groundstate is the standard 
BCS superconductor with one COMM and therefore no 
fractional flux quanta occur, (ii) For the pairing interac- 
tion (6), the groundstate for zero flux is cither a rf-wave 
BCS, or a PDW superconductor with the two COMMs 
q and — q, i.e. Q = 2, depending on the interaction 
strength V. If V exceeds a critical value V c , a first order 
transition from the BCS to a PDW state occurs. This 
phase transition originates from the competition between 
the cost in kinetic energy -Ekin and a gain in condensation 
energy E con arising from finite COMMs in the presence 
of gap nodes. The d-wave superconductor has four nodal 
points in k-space where unpaired electrons remain at the 
Fermi energy. In the PDW state, these electrons can be 
paired by A(k, ±q). The condensation energy is there- 
fore optimized by the vector q for which the nodes in 
A(k, q) and A(k, — q) are furthest apart [25]. 

In the following we investigate the flux dependence 
of the generalized PDW and its total energy E qit<l2 (p) 
for the state with the two COMMs {qi,q 2 } at T = 0. 
In particular we search for the combinations {qi, q 2 } for 
which Eq lt<ii (tp) is minimal for flux values ip = n/2, and 
we test whether combinations {qi,q2} exist for which 
Eq 1 ,q 2 { L P) is minimal at ip = n/4. The charge density is 
fixed to p = 0.8, and an additional next-nearest-neighbor 
hopping t' — —0.3t is included. For these values and 
ip = 0, a PDW state is realized for V > V c p=> 2.1 1 with 



<7i = — qi = q as the integer closest to No /A, where 
A fa 6. The value of A depends on p and typically shrinks 
with increasing interaction strength V [25]. At the flux 
value <p = 1/4 we expect a realization of COMMs with 
+ <?2)/2 = 2tp = 1/2. In particular, the combina- 
tion {q, —q + 1} is a likely candidate for the groundstate. 
However, especially for small cylinders, the favored com- 
bination of COMMs may vary with system size; therefore 
different sets {(71,92} are separately compared for each 
system size. 

Figure 2 shows the energy E q _ q+ i(l/4) as a function 
of the circumference N$ of the cylinder for V\ = 2.2 1 and, 
for comparison, E q ^ q (0) and E q _ q {\/4). The ground- 
state COMM q is the integer closest to ~ Ng/A and 
switches to the next integer when Ng increases by A. 
Since A is typically not commensurate with the lattice, 
there is a sensitive size dependence for small cylinders 
(Ng < 300) where the existence of an energy minimum 
at ip = 1/4 is not certain. For large cylinders, the total 
energy oscillates in Ng with periodic minima at Ng = nA 
and Ng = nA/2. The amplitude of the oscillations de- 
creases with 1/Ng and in the limit Ng — > 00 one finds 
that E q - q (0) = Eq t -q+x(l/4:) is always smaller than 

Eq,-q(lh). ' 

The SC groundstate around p = 1/4 is therefore re- 
alized for the combination {q, — q + 1} and E q ^ q+ i(p) 
is indeed minimal at ip = 1/4. Upon approaching 
p = 1/2, the groundstate switches to the combination 
{q + 1, — q+ 1} with an energy minimum at tp = 1 /2, and 
to {q + 1, — q + 2} around ip = 3/4 (see Fig. 3). For an 
arbitrary flux value, the combination of COMMs mini- 
mizing the energy is given by 

qi + q 2 = floor (4^3 + 1/2) , (9) 

where floor(a;) is the largest integer smaller than x. The 
new energy minima correspond to fractional flux values 
np/4 and therefore lead to flux quantization in units of 
$o/2. This result remains valid also for thick- walled, ir- 
regular or disordered cylinders, which can be modeled 
by averaging the total energy over various channels with 
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FIG. 3: Groundstate energy E(cp) = min,j li92 E qizq2 (ip) for 
the same model parameters as in Fig. 2 with Ng = 796. The 
numbers in the shaded areas indicate the COMMs qi and 
q2 in the corresponding flux regimes. The energy difference 
between the blue and the pink minima corresponds to the 
difference between the blue and pink dots at Ng — 796 in 
Fig. 2. Wether the fluxless state is formed by two odd or two 
even COMMs depends on the system size Ng. 

different circumferences, fn the thermodynamic limit, 
the energy minima at flux values $ = and $ = 3>o/2 
become degenerate, and consequently, the flux periodic- 
ity of thermodynamic properties like the supercurrent is 
<i>o/2 as well. This situation is illustrated in Fig. 3 for 
N g = 796, for which E q _ q (0) and E q _ q+1 (l/4) are al- 
most equal. The groundstate energy as a function of ip 
forms a series of parabolae with minima at ip = n/4. The 
supercurrent is obtained as 
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(10) 



and vanishes whenever the groundstate energy has a min- 
imum. Therefore, the formation of a SC state with half- 
integer valued phase winding number (gi+fft) /2 compen- 
sates the current induced by the half-integer magnetic 
flux $ = $o(9i + 92) /2 and thereby minimizes the total 
energy of the system. 

In finite systems, SC states with even and odd phase- 
winding numbers are distinct and their energy minima 
at the corresponding flux values differ [31]. The flux pe- 
riodicity is therefore strictly /ic/e, although quantization 
occurs in units of hc/Ae. We thus identify three dis- 
tinct classes of pair states on a flux threaded cylinder 
(see Fig. 3): superpositions of two odd (purple) or two 
even (red) COMMs with energy minima at integer multi- 
ples of <&()> and even-odd (blue) combinations with energy 
minima at half-integer multiples of $ - These pair states 
correspond to integer or half-integer phase winding num- 
bers of the SC order parameter. In principle, states with 
more than two different COMMs are possible as well. 
Such states will allow for fractional flux quanta of mag- 
nitude $o/Q; but they are not realized as groundstates 
in the parameter range of our model Hamiltonian. In al- 
ternative Hamiltonians the prerequisite will remain the 



existence of an instability towards nodal superconductiv- 
ity, which allows for fractional flux quantization due to 
coexisting COMMs. 

In a real flux threaded SC ring, flux always penetrates 
into the superconductor itself. The system is then no 
longer periodic in $ and the spatial variations of the mag- 
netic field induce additional COMMs. Therefore, a uni- 
directional PDW state is expected to break down even 
in weak magnetic fields, leaving a standard d-wave su- 
perconducting groundstate. Above this threshold field, 
$o/2 flux-quantization will then disappear abruptly, re- 
turning to the flux quantum $o- As described above (c.f. 
Fig. 1), in YBCO SQUIDs the periodicity of the super- 
current changes indeed abruptly at a threshold magnetic 
field [12]. The low- field $o/2 oscillations were also re- 
ported to vanish close to T c [12]. While the solutions of 
our model do not change qualitatively with temperature 
far below T c , the PDW is also replaced by a d-wave BCS 
state upon approaching T c . 

Whether a PDW state occurs in cuprate supercon- 
ductors is unsettled. However, uniaxial charge and spin 
modulations in coexistence with superconductivity were 
verified for 214 cuprates [32-35], which inevitably imply 
spatial modulations in the pair density [36-40]. In the 
typical rectangular geometry of the SQUID experiments, 
the orientation of an anticipated PDW would differ from 
the cylinder geometry discussed above. Although the 
emergence of SC states with half-integer phase-winding 
numbers and $0/2 periodicity will inevitably occur, too, 
the explicit solution of a model Hamiltonian in a SQUID 
geometry is more involved — especially for a finite mis- 
match angle between the lattices on both sides of the 
Josephson junctions. A possible analysis of this geome- 
try will require numerical simulations. 
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